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The quantum Langevin equation is derived from the Feynman- Vernon forward-backward path 
integral for a density matrix of a quantum system in a thermal oscillator bath. We exhibit the 
mechanism by which the classical, c-valued noise in the Feynman- Vernon theory turns into an 
operator-valued quantum noise. The quantum noise fulfils a characteristic commutation relation 
which ensures the unitarity of the time evolution in the quantum Langevin equation. 



. I. INTRODUCTION 

o\ : 

0^ ■ The quantum Langevin equation 

5— ( ■ mx t +jx t + V'(x t )=fj t (1) 



(N 



describes successfully the temporal behavior of a quantum mechanical point-like particle with the action 



t 



m 



A S [x]= I dt r ^ ± 2_y {xt) (2) 



> 

in the presence of dissipation The time argument is indicated by a subscript, for notational brevity. The 

dissipation is accounted for in Eq. (|l]) by the two phenomenological terms: the friction term jxt> and the quantum 
noise operator f] t which satisfies the commutation rules 

; //,./„•] = 2ii n mi ■ /') on 

■ and has the expectation values 

^g; ^{[mM+)f) = K T {t,t') , (4) 

a ■ 



with 



K T {t,t') = — [ dwwcoth-^ cosw(t-t'). (5) 

2kT 



7T 







X 

H 

The prefactor 7 is once more the friction constant, in accordance with the fluctuation-dissipation theorem. As usual, 
T denotes the temperature and k the Boltzmann constant. Remarkably, the temporal evolution governed by the 
quantum Langevin equation preserves the canonical commutation relation, i.e., quantizing [i t ,p t ] = ih with p t = mx t 
at one time, this remains true at all times 

The quantum Langevin equation (|l|) has been derived Q-^] from the Heisenberg equations of motion of the quantum 
system coupled to a thermal environment consisting of a bath of infinitely many harmonic oscillators. The coupling and 
the spectral distribution of the oscillators are chosen in such a way that the Heisenberg equations for the environment 
coordinates can be solved explicitly, and that the standard friction term arises in ([!]). The general procedure is 
described in (t^. 

The quantum Langevin equation can be used to derive Kubo 's stochastic Liouville equation JsJ for the temporal 
evolution of the density matrix of the system || . Solutions of this equation are in agreement with experimental data. 
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The same Kubo equation can be obtained |10[ from a completely different description of the quantum system due 
to Feynman and Vernon jl^,[ll| . Here the density matrix of the quantum system in an oscillator bath is represented 
by a path integral which contains a fluctuating noise variable rjt- In contrast to the quantum Langevin equation, the 
fluctuating noise in this formulation is classical. 

In ordinary quantum mechanical path integrals, there exists a simple way of going from fluctuating classical variables 
to Heisenberg operators. For the classical fluctuating noise variable in the Feynman- Vernon path integral, on the other 
hand, the relation to the quantum noise operator r)t is unknown. The purpose of this paper is to exhibit this relation 
and to derive the quantum Langevin equation from the Feynman- Vernon path integral. 



II. THE FORWARD-BACKWARD PATH INTEGRAL FOR THE DENSITY MATRIX OF A SYSTEM IN 

A THERMAL BATH 

Consider the quantum mechanical system described by the action (^) in contact with an oscillator bath. The density 
operator p t of the total system obeys the Neumann equation 

ihdtpt = [H, p t ] , (6) 

where the Hamilton operator is the sum of three terms 

H = H S + H B +H 1 . (7) 

The first term H s describes the quantum system by itself, H B the oscillators in the heat bath, and H 1 the interaction 
of the system with the bath. A formal solution to Eq. reads 

/ iHt\ „ / iHt\ 
Pt = eXP {-—) PoeXP {—) ■ (8) 

Let X denote collectively an infinite set of oscillator coordinates in the heat bath. Then a matrix element of pt has 
the time evolution 

(xX\p t \x'X') = J dxdx 1 dXdX\xX\e~ lBt/n \xX){xX\p n \x' X'){x' X'\e 1 ^ (9) 

By time-slicing the matrix elements of the two evolution operators on the right-hand side (one of which working 
forward, the other backward in time), they can be represented by corresponding path integrals H-ll^], 

with the boundary conditions xq — x, xt — x and Xq = X, X t = X, and 

(x'X'\e lBt/h \x'X') = f -Dx'T>X'e- iA l x '< x 'V n , (II) 



with the boundary conditions and X' Q = X', X[ = X'. The functional 

A[x,X] ee A s [x] +A B [X] +A I [x,X] (12) 

is the classical action associated with the Hamiltonian (^) of the total system. 

We are interested in the temporal behavior of the system coordinate Xt regardless of the behavior of the environment 
variables X t , i.e., we look for the so-called reduced description of the total quantum system. This is obtained by tracing 
out the environment coordinates in (|9|). The system and the bath may be assumed to be decoupled at some initial 
time, say t = 0. At that time, po = PoP B , and we can rewrite (^) as 

pf(x,x') = J dX(xX\p t \x'X) = J dx\dx'-JJt(xx' , x\x\)p^ {x\, x\) , (13) 

where Ut(xx > , Xxx'i) has the forward-backwardpath integral representation (the hyphen is pronounced minus) 
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U t (xx',xix[) = / Vx + Vx~ F[x + ,x-]exp ^ {A s [x + } ~ A s [x-]} 



(14) 



with x+(0) = Xi,x+(t) — x,X—(0) — x'i,x—(t) = x' . For variables with subscripts ± we place the time arguments in 
parentheses, to avoid a pileup of subscripts. The functional 

(15) 
(16) 



F[x+,X-] = I dXdX'p B (X,X') 

dX [ VX + VX_exp^ {A I [x+,X+} -„4 J [x_,X_] +A B [X+] - A B [X_}} 



0-|l3|. The forward-backward path integral over the bath coordinates is taken 

jp=x,x+( 



is Feynman's influence functional 

with the boundary conditions X + (6J~=~X, X + (t) = X,X-(0) = X' and X-(t) = X. 

The simplest bath action containing oscillators of frequency lo with coordinates X u (t) is 



A B [X] = \\di!\dw 



Xl - u?Xl 



o o 

Assuming an initial equilibrium density operator for the bath, 

p B =cM-H B /kT), 
the interaction leading to the correct quantum Langevin equation is [|3|.[13[|1 



T OO 

A^XfX] = (27/tt) 1 / 2 J dt'x J dwX^ 



(17) 



(18) 



(19) 



Since all integrals are Gaussian, the influence functional can be calculated explicitly 1 1 1 - 13 1 . The result is 



F[x-l,X-] = F o exp I — / dt'yx | exp 
\ hJ o J 

= GXP ?/ dt ' y± ) Ffl[y] ' 



2h z 



dt'dt"y v K T {t',t")y v , (20) 

(21) 

(22) 

The constant Fo = F[0, 0] in j2Cj ) is the initial value of the influence functional, where the system and bath are 
decoupled, i.e., Fo = Trp B = J dXp B {X, X); the exponential in ( ^l| ) is the dissipation part of the influence functional, 
while the functional Fn[y] describes thermal and quantum fluctuations. The two factors contribute in completely 
different ways to the Langevin equation. 

Introducing a new function pf (x, y) = pf(x + y/2, x — y/2) and substituting ( pp| ) into ( p^ ) and (|l3|), we obtain the 
evolution equation 



where we have introduced new variables 

x = (x + + x_)/2, y = x + — X-. 



Pt{ x iV)= I dx'dy'U?(xy,x'y')po(x',y') , 



where Uf (xy,x'y') is given by the functional integral 



U t s (xy,x'y') = F I VxVy e X piA T [x,y]/h 



(23) 



(24) 



with the boundary conditions xq — x',Xt — x,yo = y',yt = U- The temperature-dependent action is 
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A T [x,y] = J dt' 



% * 

myx - jyx - V(x + y/2) + V(x - y/2) + T^V K Ty 



(25) 



For brevity, we have written Kryv for the integral f dt"KT(t',t")yt» in the last term. Due to this term, the 
temperature-dependent action is nonlocal in time. Note that at t — 0, Uq (xy, x'y 1 ) reduces to Fq5(x — x')S(y — y'). 

The temporal nonlocality of the action ( p5| ) can be removed by means of a random noise variable with the correlation 
function 

(VtVt')r, = K T (t,t'). (26) 

We simply make use of the Gaussian identity 

cxp I J dt 'yKTy\ = Jvrjexp J dt' ^-^nK T 1 n + ^qy^j 

V o / o 

t 

= (exp l - J dt'y^r, . (27) 



associated with (f26|). Thus, if we agree to average all equations at the end with respect to the noise variable, we may 
replace the effective action in the evolution operator (|24|) by 

t 

A n [x,y] = J dt' [myx - jyx -V(x + y/2) + V(x- y/2) + yrf\ . (28) 
o 

In the limit of high temperatures, the noise fluctuations become large and local, 

(VtVt')n = K Ht, - 2fcT 7 ^ - t') + 0(l/T) . (29) 



Equation (|27|) shows that the average size of the fluctuating variable y goes to zero like h/^/kT. Hence the action 
( p5| ) becomes local, the potential difference V(x + y/2) — V(x — y/2) can be approximated by V'(x)y, and ( |28| ) turns 
into 



An integration by parts has been performed in the first term of the integrand in (|28|), absorbing a boundary term 
in the prefactor of the path integral (Q). We can now perform the functional integral over y and find the classical 
version of the Langevin equation ([!]) : 

mx t + jx t + V'(x t ) = T) t , (31) 

which describes classical Brownian motion. The classical behavior at large T is a consequence of ft, and T appearing 
in the combination fttu/2kT in (|^), so that T —> oo is equivalent to ft —> 0. 

There have been attempts to include quantum effects into this classical equation by simply replacing the local 
noise correlation (p9h by the finite temperature one (p^), while maintaining the linear appr oximation V'{x)y to 



V[x + y/2) — V(x — y/2). The result has been called a quasiclassical Lanyevin equation [|14 15|. Obviously, such an 
approximation can be reasonable only for nearly harmonic potentials (see also the remarks in Ref. [ [To| , p. 589). 

It is surprising that by converting the variables Xt and r/t in (|3l]) into operators in ([l]), all terms in the expansion 
of V(x + y/2) — V(x — y/2) in powers of y can be accounted for, as we shall prove. 

III. KUBO'S STOCHASTIC LIOUVILLE EQUATION 

Equations ( |23| ) and (|24|) resemble the evolution equation of a Schrodinger wave function whose role is now played by 
the density function pfjx,y). This analogy is helpful in deriving a Schrodinger-like differential equation for pf . This 
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equation is most easily found by going over from a Lagrange-type path integral (|2J) to a Hamiltonian one, expressed 
in terms the position and canonical momentum variables fllql . The effective action goes over into 

t 

A T [x, y] -> A T [p x ,Py, x,y} = J d£ (p x x + p y y - H T ) , (32) 



where 



H T = - {p x + jy) P v + V{x + y/2) - V(x - y/2) - ^-yK T y (33) 
m 2n 



plays the role of a Hamiltonian. For the noise-dependent effective action (f28| ) we find an analogous action 
A v [pxiPy,x,y\ involving the Hamiltonian 

H v = -(p x + 7 y) P V + V(x + y/2) - V(x ~ y/2) - yr,. (34) 
m 



Using the effective action (p2j), we obtain the alternative path integral representation for the time evolution operator 



Uf(xy,x'y') = F J Vp x T>p y Vx'DyexpiAT\p x ,Py,x,y\/Ti . (35) 

A similar equation holds for a fixed noise with the action Ari\p x ,Py,x,y]. 

Since the action A v [p x ,p y , x, y] is local in time, we conclude that for a fixed noise rj, the associated noise-dependent 
{noisy) density matrix pf[rj\ obeys the Schrodinger-like equation 

thpfM^H.pfin] , (36) 

with H v being the operator arising from (|34|) by substituting p x —> —ihd x , p y — > —ihd y . The noise average of the 
solution to ( |36"| ) is the density matrix of the system: 

~p s t {x,y) = {p s t [ri\) v . (37) 

Equation (|36| ) supplemented by (^6|) and (^) is called Kubo's stochastic Liouville equation 0], which we have thus 
derived from the Feynman- Vernon forward-backward path integral. 

Note that there is no analogous procedure to obtain a differential equation for the density matrix pf(x, y) from the 
Hamiltonian (|33| ), due to the nonlocality of the i^r-term. Only in the limit of large temperatures, when Ht becomes 
local due to (|29|), there exists a Schrodinger-like equation which is the Fokker-Planck equation at finite friction. 

It must be pointed out that the transition from the Hamiltonians Ht and H v to their operators has an ordering 
ambiguity in the term r yyp y /m in Ht and H n |l6| . In the time-sliced path integral we must decide whether to 
write jyuPyn/m or 72/ n p yn _i/m where the subscript n numbers the time slice. Since the ordering is independent of 
temperature, we resolve the ambiguity in the limit of large temperatures: from the well-known Fokker-Planck equation 
we determine the correct operator ordering of the term ^ypy/m to be —ih^ydy/m. 

IV. CLASSICAL NOISE VERSUS QUANTUM NOISE 

The noise in Kubo's equation is a c-number, and its relation to an operator-valued noise of the quantum Langevin 
equation (Jl]) has been an outstanding puzzle, as pointed out in the Introduction. To find this relation, we observe 
that it is possible to remove the temporal nonlocality in the initial effective action ( p5| ) just as easily with the help 
of an operator-valued noise r). This was, in fact, how the Feynman- Vernon path integral was derived in Ref. [fj"3"[ [see 
Eq. (18.162)]. An operator- valued noise possessing the properties @-(§) can be chosen as a sum of all Heisenberg 
operators of the oscillator velocities. Explicitly: 

oo 

fjt = V^lNdt J due^ Bt / h X u e-^ Bt ' n (38) 



oo 

= *(7A) 1/2 / dcoV^ (e*-*6t - e-^aj) ; (39) 
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where and a w are time- independent creation and annihilation operators with the usual commutation rules: 

[fio,, a u ,] = o, [at , a* ,] = o, [a*, at,] = m(cc - a/) . (40) 

The noise operator satisfies the commutation rule (||). The correlation function dJ) follows if we define fj- 
averages as bath averages: 



([%,Vt'}+)f, = Tr (pH[fjt,f)t']+) 



(41) 



Thus the operator (J3q) has precisely the properties of the quantum noise variable in the quantum Langevin equation 

©• ... « 

In the operator representation of the noise variable, the influence functional is given by [|17| 



x + ] = Tr 



/5 T C exp- / dt'xcvc 
** c 



(42) 



where C is a closed-time contour encircling tightly the interval [0, t] in the complex t plane and Tq is the ordering 
operator along this contour. The subscript C distinguishes the upper and lower branches of the integration contour, 
where xc(t) is equal to x+(t) and X-(t), respectively, whereas the operator fjc is equal to fjt on both branches. 
Introducing the symbol T" 1 to denote anti-time ordering, Eq. (E2h takes the more explicit form 



F[x-, x+] = Tr 



Pa 



— J dt'x-f) J T exp I — I dt'x+f} 



(i 



(43) 



Assuming the explicit noise representation ( pq ) through the heat bath operators we observe that the matrix element 
(X\T exp(« Jq dt'x+fj/h)\X) and the corresponding matrix element of the anti-time ordered exponential in (^3| ) are 
given by the path inte gra ls over X + and X- in ( |l6| ) , respectively, whereas the trace in ( ^3| ) corresponds to the integrals 
over X, X and X' in (Q) . 

We now make the key observation that will allow us to derive the quantum Langevin equation (|l|): The fluctuation 
part of the influence functional, defined in Eq. (f2~l|), is obtained from (^3|) by setting x = 0. Then both exponentials 
in ( p3| ) carry the same argument (i/2h) J dt'yfj . When expanding the product of the two exponentials in a power 
series and performing the proper time orderings, the bath average can be rewritten as 



F R [y] = F[y/2,-y/2] = Tr 



Texp I ^ / dt'y t adfft ] 



(44) 



where ad^t is the adjoining operator associated with the noise operator r) t . This is defined as follows. When acting 
upon an arbitrary operator O, the adjoining operator is equal to half the anticommutator: 



adf, t d=[f, t ,d} + /2. 

The bath average (|44| ) corresponds to defining a noise average of an arbitrary functional ©[ad?)] by 

<C[ad fj])f ee Tr 0[adf)]p* . 



(45) 



(46) 



Note that the adjoining operator ad 77 acts also upon . The ket superscript ad emphasizes this fact which is in 
contrast to averages with respect to the ordinary noise operator defined as in (^). For example, the correlation 
function is calculated as 



(ad ?? t ad rj t >)f = ^Tr %, [vt>,Po}+}+ = ^Tr po [fjt, Vt']+ = K T (t,t') 



(47) 



The adjoining operator adr) t has an important property crucial to the further development: Although the original 
noise operator fj t has a nontrivial commutator (^) with itself at a different time, two adjoining operators ad fjt and 
ad?7 t / commute with each other, so that they can be treated as c-numbers. Indeed, for an arbitrary operator O we 
verify that 



[adr) t , ad fft']d = [vuVt'P + 6%>fk\ = 0- 



(48) 



G 



In all such calculations, we may treat the quantum noise ad fjt as if it were a classical noise rjt- 

Furthermore, it is clear from the action ([[7]) that ad fjt is a Gaussian noise variable, just as the classical rjt- The 

correlation functions of any number of operators ad fj t evaluated in the average ( |46| ) are completely specified by their 

two-point function ( [47j) following Wick's expansion rule. Since the expectation values of the products of two r\t or 

of two ad fjt are identical, < r\ t r\t' >r/=< &dfjt&dfjt> all correlation functions must be identical as well — quantum 

and classical noises are completely equivalent. 

Substituting ( py| ) into and repeating the arguments leading to Kubo's equation in Section III, we again end 

up with a modification of Kubo's stochastic Liouville equation, in which the noise is operator- valued, i.e., rjt in ( |3^ ) 

and d37| ) is replaced by ad fj t ■ 

Having shown how the operator-valued quantum noise turns into a classical noise in Kubo's stochastic Liouville, 

we are ready to derive the quantum Langevin equation (|l|). 



V. QUANTUM LANGEVIN EQUATION 

As the first step, we derive the evolution equation for the system's density matrix operator pf whose matrix 
elements are defined by dl3|). Its operator representation is pf — Tr BPt where Tr^ denotes the trace over bath 
degrees of freedom. Returning in Kubo's equation ( |36| ) to the initial forward-backword variables x± of (p2|), we find 
the evolution equation for the noisy matrix elements associated with (filf): 

iWtpf[x + ,xJ) = — Hf_ + ^-(x+ -x-)(p+ -p-) - (x+ -x_)ad fjt] pf(x+,x-) , (49) 
L 2m J 

where p± = —ihd/dx±, and H± = H s (p±, x±) is the system Hamiltonian expressed in terms of the forward-backward 
phase space variables. 

Of course, this equation can also be obtained directly from Eq. (g) . The Hamilton operators in the two exponentials 
yield directly the difference Hf — in ([l9]) , whereas the other two terms result from the dissipation and fluctuation 



parts of the influence functional in the path integral ( 14 ) 



Let \x) and \p) denote eigenvectors of the canonically conjugate system operators x and p, respectively. Then 
pf(x +> X-) =< x + \pf \x- >. We have the following identies 

(Hl-H s )(x + \p?\x_) = (x + \ [H s ,pf] ; (50) 
{x + p++x-p-){x+\pf\x-) = (x+\ (xppf -pfpx) \x-) ; (51) 
+ x + p-)(x+\pf \x-) = (x+\ (ppfx- xpfp)\x-) , (52) 

which can be proved by inserting suitable resolutions of unity, J dp\p)(p\ = J dx\x)(x\ = 1, between the operators on 
the right-hand sides. Substituting Eqs. fl50|) (^2|) into (E9h, we find the evolution equation for the operator pf [ad fj] 



ihd t pf = [H s , pf ] + [jp/m - fj t , pf ]+] . (53) 

The average ( fl6| ) of a solution to (^3|) gives the density matrix operator of the system 

pf = (pf [ad fj])^ = Tr B pf [ad fi}^ (54) 

Now we turn directly to a derivation of equations of motion for system operators. Let us recall first that in 
quantum mechanics there are two equivalent representations of the mean- value (O)t of any dynamical variable O 
(with or without explicit time dependence) at a time t 

(0) t = TrOpt = Tr6tp , (55) 

where pf is the density matrix operator at the time t, and Ot is the Heisenberg operator coinciding with O at t = 0. 
The first equality expresses the mean-value in the Schrodinger picture, while the second one determines it in the 



Heisenberg representation. Equation (55) can also be regarded as a definition of the Heisenberg operator Ot- The 
consistency of this definition is guaranteed by the time independence of the commutation relations [O, (D'\ = [Ot,0' t ] 
for any two operators. This, in turn, is a consequence of the unitarity of the temporal evolution of the density matrix 
described by the Neumann equation (^|). 

The evolution of a quantum systems with dissipation is determined by a generalization of the Neumann equation 
(l53|). We shall demonstrate in analogy with (Ej|), there exists a consistent definition of noisy Heisenberg operators, 



7 



if the noise obeys the commutation relation (||). For any product of system operators O s 
noisy Heisenberg system operators is defined by the equality 



nLi^ (/c) ' a product of 



Tr s 6 S pf[adry] 



Tr s J] G^pf [adfj] = Tr s ]J d[ k) [f}]^ = Trgdffflpo 



(56) 



fc=i 



k=l 



where the noisy density matrix operator satisfies Kubo's stochastic equation (p3|), and Trs implies a trace over system 
degrees of freedom. One can think of the relation (|56| ) as a mapping of an algebra of system operators at the time 
t = onto their algebra at a time t. Taking a trace of the evolution equation (|53"1), we see that the unit operator / 

remains invariant under the mapping (|56|), i.e., It[fj\ — 0, so that It[fj] = I. The consitency of the definition ( |56| ) is 
guaranteed by the time independence of the commutation relations 



[6 s , 6' s ] = [6f[f)],d' t s 



(57) 



for any two Of \f]\ and C^ 5 ^]. At t = 0, equation (57) is trivially fulfilled. To prove the time independence, we derive 
the noisy Heisenberg equation for Of . Differentiating (^6j) with respect to time and making use of both (|53| ) and 
([56]), we obtain 



ifi6t S 



Pf,H t s 



jpt/m-fit, 



Pf, x t ] 



(58) 



where Hf — H s (pt, Xt). For a noise operator satisfying the commutation rule (|^), this equation leads indeed to (|5 
(see Q). Equivalently, we demonstrate that Eq. ( |58| ) coincides with the Heisenberg equations of motion of the entire 
system (the particle in a heat bath), comprising the quantum Langevin equation (PP). 
For an arbitrary system operator we have 



ihd t s = [6f,H t } = [6f,H t s ] 



2 I 7T 



1/2 



dL0Xu U [6f,X t ] 



iKX wt = [X wt ,H t ] , ihP^ = [P ut ,H t ] , 
where H t is the canonical Hamilton operator corresponding to the total action ( |l2| ) defined by (||), ( |l7|) and ( |l9| ) 



H t = H? + Hf + Hi = H? + - dw 



Put - I — 
7T 



27 



1/2 



XL 



2 v2 



uj z X 



(59) 

(60) 
i: 

(61) 



The w-integral in this Hamiltonian produces a linearly divergent term oc i 2 /2, the "frequency shift term" discussed at 
length in Ref . . The divergence is due to the specific form of the interaction (|l9| ) which generates the desired time- 
independent friction coefficient 7 in the Langevin equation . The same divergence appears also in the Hamiltonian 
form of the forward-backward path integral (IXJ) . The divergence is canceled by an equal divergence in the dissipation 
part of the influence functional arising from the coupling to momenta of heat bath oscillators (the term ~ xP u in (|6l|)) 
p0| , p. 602. For this reason, Kubo's stochastic equation has only finite terms. This pleasant cancelation is absent if 
the system is coupled directly to the positions of the bath oscillators rather than to their velocities, as in (|l9|). In 
that case, the initial potential V(x) requires a divergent counter term in order to obtain a finite Kubo equation. The 
entire issue is, of course, somewhat academic since it is a consequence of insisting upon a constant friction, which is 
physical only at low frequencies, much lower than the collision rate in the system. For larger frequencies, the physical 
behavior is described by the Drude friction 7(0;) oc 1/(1 + uj2 / uj d)^ an d there is no divergence at all M, Q. 
Imagining the presence of a frequency cutoff in (|6l|), we solve the oscillator Heisenberg equations (|6G|) and find 



1/2 



— I / dujX^t = ix t - fit = ypt/m- fit 



(62) 



where we have used the representation (39) for the noise operator. Substituting 
(|58|) we conclude that 



3) into (J59|) and comparing it with 



6f[fj] = e ^i^6 s e- liit i h = Of 



(63) 
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Thus, the solutions of the evolution equation ( pq) coincide with Heisenberg system operators. Hence, the temporal 
evolution is unitary and ( ^7|) is satisfied. Note that the system Hamiltonian is no longer an integral of motion, 
dHf jdt 7^ 0, in contrast with the Hamiltonian driving the ordinary Heisenberg equations. 



Setting Of equal to it and p t in Eq. ( p8[) or, equivalently, in Eq. (59), we obtain the Heisenberg equations of 
motion: 

it = Pt/m , (64) 
Pt = -V'(xt) -jpt/m + f) t , (65) 

which are equivalent to ([!]) 

Equation ( |63| ) ensures the unitarity of the temporal evolution and the time independence of the canonical commu- 
tation relation [£t,f>t] = iti. This guarantees also that correlation functions of noisy operators within the fj average 
agree with the ordinary quantum-mechanical correlation functions of the system described by the Hamilton operator 

Tr s p§{[6? [fj], &f = Tr p [6f , 6' t f] . (66) 

The property (||) of the quantum noise is crucial for the unitarity of the time evolution of the Heisenberg operators. 
To see what happens if the noise were commutative, we replace ad fj in (|53|) by a c-valued Gaussian noise 77, which 
we are apparently allowed to do after the above observations on the equivalence of averages. Then we can again 



define the noisy Heisenberg representation by (|56| ) and derive (58) in a similar way. However, we can never prove the 
consistency of ((5^) and ([56]); Eqs. (E^) and (|63D are no longer valid. For a commutative noise, the violation of the 
canonical commutation relation as time proceeds is seen from (|64|), (|65|), most simply for V = 0. This break-down of 
the unitarity in the temporal evolution would prevent us from identifying (|5^) with the quantum Langevin equation. 

Thus we have shown that the Feynman- Vernon path integral description of quantum Brownian motion is completely 
equivalent to a description in terms of the quantum Langevin equation (|l|). Both descriptions can be used on equal 
footing to study dissipative processes in quantum mechanics. 
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